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Abstract 

We study the free energy of the Parabolic Anderson Model, a time-continuous model 
of directed polymers in random environment. We prove that in dimension 1 and 2, 
the free energy is always negative, meaning that very strong disorder always holds. 
The result for discrete polymers in dimension two, as well as better bounds on the free 
energy on dimension 1, were first obtained by Hubert Lacoin in IITOl , and the goal of 
this paper is to adapt his proof to the Anderson Parabolic Model. 

Keywords : Parabolic Anderson model; Brownian Directed Polymer; Free energy; 
Strong disorder; Coarse graining. 

1 Model and known results 
1.1 The Parabolic Anderson Model 

The general model we consider in this paper is defined in terms of a random trajectory 
in a random environment : 

• The random trajectory : let {(Xt)t>o, -P") be a nearest-neighbour random walk on 
the lattice Z"' with rate of jump k starting from 0. Specifically, we let (M, 5") be the 
path space of cadlag trajectories from to Z^, and the process {Xt, P"^) has for 
infinitesimal generator kA'^, where A"' is the discrete Laplacian. Sometimes we 
will need to start the random walk from a point other than 0, let {{Xt), P^) be the 
walk starting from the point x. 

• The random environment : let {{B^{uj))t>o,x E Z'^) be a family of independent 
Brownian motions defined on a probability space (fi, 9, Q), one for each vertex of 
Z*^. For t > 0, it is natural and convenient to introduce the sub-cr-field 

5t = a{B-:,se [0,t],xGZ'^}. 
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• The Gibbs measure : for a time horizon t we define the following Hamiltonian : 

For any t > 0, we define a probability measure /ij' '^ on the path space J) 

where /3 > is the inverse temperature and 

Z^^^ = P^ [exp(/3/7,(X)-t/3V2))]. 
Since the polymer measure is parametrized by the environment u, it is random. 

This model has been studied in 141 with the point of view of intermittency, and in 
||9| when interacting with a particle system. It has also been studied as an example of 
polymer model in |fT|, fl\. In [13J Moriarty and O'Connell studied a perfectly assymet- 
rical version of the 1-dimensional model in which they could prove exact formulas. 

In this paper, we investigate the asymptotic behavior of the partition function ''^ 
when t — oo. 

Let us finish the definition of the model with some remarks on scales. The process 
(Xt, P'^) has the same law than (X^t, P^), and if Bt is a Brownian motion, fi^t has the 

1,— 

same law than ^/kBi. Therefore, Z^'^ has the same law than Z|^^^ , and we only need 
to study the case of k = 1, and the rest is obtain by rescaling. 

We will denote by P, /if, Zf,... the quantities P'^, li^'^ , Z^'^ with k = 1. In most 
case, when the value of (3 is not ambiguous, we will only write ji-t and Zt. 



1.2 The partition function 

Let us fix a /3 > and begin by stating that the partition function Z^ is a martingale on 
(M, S, Q). For any fixed path X, the process {Ht{X)}t>Q is himself a Brownian motion, 
and {exp(/9ift — t(5'^ /2)}t>Q is its exponential martingale. Therefore, the partition func- 
tion Zt is a mean-one, continuous, positive martingale on (M, S, Q) with respect to the 
filtration (St)t>o- In particular, the following limit exists Q-a.s. : 

Zoo = lim Zt 

t—^oo 

Since exp(/3i7j) > 0, the event {Z^o = 0} is measurable with respect to the tail cr-field 

t>i 

and therefore by Kolmogorov's — 1 law, 

Q{Z^ = 0} = or 1. 

In the case Z^o > 0, we say that weak disorder holds, and in the case Zoo = that strong 
disorder holds. 
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1.3 The phase transition 



Theorem 1. There exists a critical value (3c = (3c{d) e [0, oo) such that : 

• if /3 < /3c then > a.s. 

• if /3 > /3c then Z^o = a.s. 

Moreover, ifd>3 then /3c{d) > 0. 

A lot of work has been made to study these two different phases for random walk 
in random environment models, and it has strong links with the overlap of two inde- 
pendent trajectories in the same environment to. If we call Jt the expected overlap of 
two independent replicas : 



For precise results on the localization of the polymer, favorite end-point and favorite 
trajectory, I invite you to consult [2] and f5\. 

Empirically, in the strong disoder case, the polymer will be compelled to follow a 
most attractive path, while in the weak disorder phase, the probability will be more 
fairly distributed amongst all possible trajectories. This is why this phase transition is 
sometimes called the transition from delocalized phase to localized phase. 

It is also interesting to know when very strong disorder holds, meaning that Zt con- 
verges to exponentially fast. It is a well-known fact in the PAM that the limit 



exists a.s. and in for p e [I, oo). It is called the free energy of the model. For more 
results on the free energy, consult [8], and read [7J for the large deviations point of 
view. 

1.4 Main theorem 

In this paper we prove that very strong disorder always holds in dimension 1 and 2. 

Theorem 2. In the Parabolic Anderson Model for dimension 1 and 2, for all /3 > 0, there exists 
c > such that 




then, roughly speaking : 










= lim -\nZ, 



Zt = 0{e-^') g-a.s. 
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Moreover, we have the following upper bounds for small (3 (ci, C2 are positive constants): 



In other terms, as t — )■ oo, the partition function grows exponentially slower than 
its expectation. 

Remark 1 : In higher dimensions, as a corollary of theorem 1, the free energy ^(1, /3) 
is constant equal to on the (non trivial) interval [0, l3c{d)) (at least). 

Remark 2 : A lot of work has also be done to study the asymptotics oi ip{l, /3) for (3 
big. It is well known in the Parabolic Anderson model (see HJ] |3l |5l ||8l) that when 

/3 cx) 



where a is a constant depending on the dimension. 

Remark 3 : In 1131 , O'Connell and Moriarty studied a perfectly assymetrical version 
of the 1-dimensional model : the random walk Xt always jumps in the same direction. 
They were able, thanks to reversibility properties of this new model, to prove an exact 
formula for the free energy. The asymptotics of this free energy when (3-^0 are — §4 + 
0(/3^), which is the order found in our theorem. Moreno in [12] also showed that an 
assymetric Parabolic Anderson Model with very high drift had the same behavior than 
O'Connell and Moriarty's perfectly assymetric model. 

The goal of this paper is to adapt the proof of the same results for directed polymers 
in random environment for discrete times, which can be found in [ilOJ to the Parabolic 
Anderson Model. 

The same proof can be adapted for other directed polymer models : a Brownian 
directed polymer in a Poisson environment as seen in 161, the Linear Stochastic Evo- 
lutions model introduced by Yoshida in |14l[T5| . Similar techniques are also used for 
non-directed models, Lacoin shows in |fTT| that the number of self avoiding walks of 
length n in the infinite cluster of the supercritical percolation in is exponentially 
smaller than its expectation. 

2 Proof in dimension 1 

Before starting, we sketch the proof and how it is decomposed in different steps. These 
steps are the same as in the section 3 of IITOl : 

1. We reduce the problem to the exponential decay of Q[Zf] for some 9 e (0, 1). 

2. We use a decomposition of Zt by splitting it into different contributions that cor- 
responds to trajectories that stay in a large corridor. 

3. To estimate the fractional moment terms appearing in the decomposition, we will 
change the measure Q around the path. 




for d = 1 
for d = 2 
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4. We use some basic properties of a random walk in Z to compute the expectation 
with the new measure. 

2.1 The fractional moment method 

To show the exponential decay of Zt, we will use the trick known as the fractional 
moment. We want to show that 

1 

lim -g[logZt] < -c < 

t— ^oo % 

but it is not easy to handle the expectation of a log. So we pick Q e (0,1), and by Jensen 
inequality : 

Q[log = ig[log Zl\ < i log Q[Zf] (1) 

We are left with showing that Q[Z^] = 0(6""'*), which is easier. In the rest of the proof, 
we will take = ^. 

2.2 Decomposition of Zt 

We fix n to be a very large number. The meaning of very large will be seen at the end 
of the proof, . The number n will be used in the sequel as a scaling factor. We also take 
n such that y/n is an integer. In the following, when we define n, the notation n = f{/3) 
will mean that n is the smallest squared integer bigger than If > 5, then 

m <n< 2/(/3). 

For z e Z, we define = [z^/n, {z + l)^/n) so that the are disjoint and cover Z. 
For T = ran, we sort the paths {Xs)se[i:i,T\ of the random walk by the position of the 
points (Xm)™ 1. More precisely, we look at the m-uplet {zi, . . . , Zm) such that Xin G hi 
for i = 1, . . . , m. 

Now we can decompose Zt : 

^ Z{z^,...,Zm) (2) 

{zi,...,Zm)&'^ 

where 

In the sequel, for shorter notations, we will write Z— {zi,...,Zm) and use the con- 
vention that zq = 0. 
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a trajectory for ^-1,2,0,3,0 



Now we use our 6* and the inequality ^CZli^if ^ '^i (which holds for any count- 
able collection of positive real numbers) to get : 



2.3 The change of measure 

We fix one Z= (21, ... , z^). We define m blocs {Jkf^zl and their union : 

Jk '■= {(s, y) G M+ X Z^, kn < s < {k + l)n and |?/ — 2;fc_i-\/ri| < C'lV^} 
where Ci will be a big constant. 



(3) 



(4) 




the boxes Jk 

Let 5 > be a small real number (we will fix its value in ((6])), we define a new 
measure on the environment : 



dQz 
dQ 



m-l / /"^ \ 



By Cameron-Martin-Grisanov theorem, under this new measure, the processes {B^).j.^-i 
are still independent and verify the following equations : 

where the {Bf)xez are independent Brownian Motions under Qz- 
We use Holder's inequality to get an upper bound on Q[Z2] '■ 



Q[Zz] = Qz 
< Qz 



dQz 
dQ 



dQ 



-1 



Qz [Zz]' 



(5) 



we can compute the first term of (O, it is standard Gaussian calculus 

-0 -1 i-e 



Q2 



dQz 
dQ 



Q 



dQz 
dQ 



exp 
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We know that the size of Jz is approximately 2Ciny/nm. We will fix the value of 5 such 
that this quantity is small (recall that 9 = 1/2) : 



dQz 
dQ 



let 6 = C;^^\-^/\ then Qz 

2.4 Upper bound on the second term 

Using dll) with ^ and ©, we obtain 



1-6 



(6) 



(7) 



Under the measure Qz, all the brownian motion in the boxes Jk have a negative drift 
—5, so if Xt stays in the boxes, Ht{X) is a gaussian variable of mean —T5 and variance 
T, and 

Qz[exp(/3ifT(X) -T/3V2)] = e"^^^. 

Now recall the value of 5 from we can take n big enough such that this quantity is 
really small (C2 will be a big constant): 



let n = ClC2p-\ then e"^^" = (e-'^^^y/"' < e'^^^' is really small 



(8) 
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Now we have to consider that the trajectory may exit the box at some time. Let 
e > 0. We note Tj^ = Tj^{X) the time spent in Jz by the trajectory X, and Tj^ the time 
spent in the box J^. 

Qz\Zz] = QzP[expi(3HTiX)-Tpy2)l{x^^^j^^,=i,...,m}] 
= P[exp(-/35TjJl|x„e/.,,i=i,...,m}] 

"m— 1 

n exp(-/35TjJl|x,„G/.,,i=i,...,m} 



p 



k=0 



A trajectory verifying Xi„ G J^. for i = 1, . . . ,m can be seen as m slices of trajectory, 
respectively starting from one point of J^.^^ and finishing at a point of J^-, which are 
similar to slices starting from somewhere in Jq and finishing in J^- So we can 

bound the contribution in the above expectation by maximizing over the starting point 
X in Jo, and we get the following upper bound 

1 

Qz[Zz]< Y[raaxP.^ exp(-/35TjJl|x„e/.,-.,_ J 

k=0 ^ ° 

where is the probability under which the trajectory starts from x. 
Now that each term of the sum ((7| is split in m independent terms : 

»n— 1 

zez™ zez™ A:=o ° 

VmaxP, [exp(-/35TjJl{x„e/4]' 

Basic properties of the random walk on Z tell us that the probability G Iz] 

decreases like exp(— when z — )■ ±oo, so we can find a R big enough such that 

J2 maxP, [exp(-/35TjJl|x„6/4]' < ^- (9) 

\z\>R ^ ° 

Now we use the following trivial bound for the remainig terms : 

V max P^ [exp ( -/35TjJ 1 e/.}] ^ 

|2|<iJ 

< {2R + 1) maxP, [exp(-/35TjJ]' . 

We now estimate roughly the time Tj^ spent in the box : Tj^ = n ii the trajectory 
never leaves the box, and is bounded from below by if it does : 

Px [exp{—(]5Tjg)] < exp{—l36n) + Pr[the trajectory leaves the box] 
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max \Xt\ > (Ci — l)\/n 

te[0,ra] 



Then we take Ci big enough such that 
P 

and we just have to inject © and ((T0|) in (O and we get : 



and taking C2 big enough, this is smaller than e To conclude, we put this result in 
(d) and © : 

3 Proof in dimension 2 

Once again we sketch the proof before starting : 

1. We reduce the problem to the exponential decay of Q[Zf] for some 6 e (0, 1). 

2. We use a decomposition of Zt by splitting it into different contributions that cor- 
responds to trajectories that stays in a large corridor. 

3. To estimate the fractional moment terms appearing in the decomposition, we 
will add a term that will handicap the environments that have a lot of positive 
correlations around the corridors corresponding to each contribution. We define 
this handicapping term in such a way that the new measure is not very different 
from the original one. 

4. We use some basic properties of a Random Walk in to compute the expectation 
with the new term. 

The steps 1 and 2 are exactly the same than in the 1-dimensional case, refer to the 
sections 2.1 and 2.2. We define the interval as follows : 

z = (a, b) e 1?, then = [oA/n, (a + \)\fn) x [fo^n, (6 + V)\fn). 

We obtain the following bound : 
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3.1 The tweaking of the measure 

We jfix one Z= {zi, . . . , z^)- We define m blocs {Jk)k=o and their union Jz : 

Jk := {{s, y) eR+x 1? , kn < s < {k + l)n and \y - Zk-iVn\ < C3^/n} (12) 



where C3 will be a big constant. In the following, I will abuse this notation and say that 

ye Jkif\y — Zk\fn\ < Cj,^ {ie y is in the vertical component of Jk). 

First of all, let me go on a little tangent. If B and 13 are two brownian motions, let 
us define what I will call the correlation between B and B in [0, T] pt{B, B) : 

Pt{B, B) = (^j^ dB^j dBs + I dB^j dB,. 

This quantity is easier to understand if we consider that the Brownian motions are 
constituted of a sum of tiny discrete increments A„Ss = ^^5+ 1 — Bg. Then the quantity 

n 

becomes 

Pt{B, B) = lim V A„S, A„ B^. 

n— >-oo ' ^ 

< s,u < T, s ^ u 

' n 

So if the quantity pt{B,B) is big, it means that there are a lot of discrete increments 

{AnBg, AnBu) that are the same sign, or in rough term that B and B are both "mostly 
increasing" or both "mostly decreasing" (actually, it is easy to compute that pt{B, B) — 
Bt ■ Bt, but I wanted to stress the "tiny increments" vision). 

The quantity that will measure the correlations of the environment in one box Jk is 
a more complicated version of what we just did : 

Rk = Rk{u^) =J2J2 / ns,=c)iu,y)dBl dB^, (13) 



where 



Vis,.)iu,y) . iQQCsC,nV^{\s -u\ + iy ^ ^ 

The quantity Rk has mean zero, and, roughly speaking, it is positive when there is 
a lot of close brownian motions which are "mostly going the same way". To compute 
the mean and variance of Rk we simply use that the B!^ are independent Brownian 
motions, on the same filtration St/ and then 



Q 



J f{s)dB^j U g{u)dBlj = J f{s)g{s)ds if x = y and otherwise. 



p{k+l)n PS 

VavQiRk) = E E / / ^(i 
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Let us compute this term : 

VavQiRo) -EE// 

^^T„,.^jJ0 Jo 



1 



{\x-y\<C4y/\s-u\} 



xeJo y&Jo 

en PS 



xeJo 



WCiCjn^ logn(|s - u\ + 1)2 

-duds 



duds 



4^2 


s — u 




W^ClClnHognils 


— u 


+ 1)2 



4^2 


s — u 




104C32C42n2 1ogn(|s 


— u 


+ 1)2 



duds 



AC! 



< 



WCiCjn^ logn 



lO'^ClClnHogn 



WClClnHogn 
< 1 



\og{n)ds 



Let Khea large constant. One defines the function fx on IR to be : 
We define gz function of the environment as 



gz{uj) := exp ^E (^t) j 



So if in a box the quantity is too big, we affect a weight e ^. 
We use Holder's inequality to get an upper bound on Q[Z2} : 

Q[zl] = Q[gz{uj)-<^{gz{co)Zzy] 

< Q[gz{ujr^Y-^Q[gz{u)Zz]' 

The block structure of gz allows to express the first term as a power of m. 

e 



Q 



e 

' i-e 



Q 



exp 



i-e 



h (Ro) 



We have already seen that Q[Ro\ = and VarQ{Ro) < 1. In consequence, 

Q [Ro > expK^] < exp{-2K^), 

hence 



Q 9z{^) 1- 



e 



< T 



< 1^1 + exp ( -K - 2K' 

if K is large enough. We are left with estimating the second term 

Q [gz{uj)Zz] = QP [gz{u^) eM^Hr - r^V2)l{x,„e7., , i=i...m}] 
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For a fixed trajectory of the random walk X, we consider Qx the modified measure 
on the environment with density 

^:=exp(/3i7T(X)-T/3V2). (19) 

Under the measure Qx, the -Bf are still independent processes on the filtration S^, 
and 

= + Pl{Xt=x}dt, 

where the i?f are independent brownian motions under Qx- 
With the measure Qx, (HH) becomes 



Q[gz{uj)Zz] = PQx[gz{uj)l{x.„^i^^, 



i=l...m\ 



Now we try to bound this term by something that could be decomposed on the blocks 

Jk- 

As we did in dimension 1, we slice the trajectory X[o,t] in m bits of trajectories start- 
ing from somewhere in Jq and finishing in J^^.^. j. So we can bound the contribution 
in the above expectation by maximizing over the starting point x in Jq, and we get the 
following upper bound 



m 



Q [gz{uj)Zz] < n max P.Qx exp (/^ (i?o)) j ■ (20) 

i=l 

Using this with ((17)), ((TSl) and ([TT]) we get the inequality 

P[Zf] < 2"^(i-'^) ( maxP.Q. [exp (Po)) l{x„e/4]' 
Therefore to prove the exponential decay of (^[IVf ] it is sufficient to prove that 

5^ S^""^^ ^^^^ ^^^^^ l{^ne/.}]'' 

is small. 



3.2 Computation of the second term 

We fix some £ > 0. Asymptotic properties of the random walk guarantee that the 
probability Q[Xn G /^] decays like e"'^'^/^ when \z\ — >■ oo. So we can find R such that 

V max P^Qx [exp Uk {Rq)) l{x^eh}] ' < V max P,[X„ G /J'' < e (21) 

\z\>R \z\>R 
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To estimate the rest of the sum, we have the trivial bound 



E 

\z\<R 



maxP^Qx [exp {fx (Ro))] l{x„e/.,-., J 



< 4i?2 max P,.gx [exp {/k (i?o))] 



(22) 



(23) 



Therefore it is sufficient for our purpose to check 

max P-,Qx [exp {fx (Ro))] < ^' 
xeio 

for some small e'. 

Now we have to find an upper bound for 

P.Qx [expifKiRo))] 

that is uniform in x G Jq. To study the left-hand term of the equation ^TT} , we showed 
that under Q, Rq had a mean and a small variance, so fK{Ro) was with very high 
Q-probability Now we will show that under Qx, Ro has a really high mean and a small 
variance, so /a'(-Ro) = —K with high Qx-probability 

We fix 5 > small. First, we force the random walk to stay in the zone where the 
environment is modified by writing 

max P^Qx [exp (fx (Ro))] < 
xe/o 



Po 



max \Xs\ > (C3 — l)\/n 

s£[0,n] 



+ maxP^Qx exp (/^^ (i?o)) l{Xro„,cjo} 



(24) 



and we can take C3 big enough such that Po[max[o,n] > (C'a — l)y/n] < S, so we can 
concentrate on the trajectories that are inside the box Jo- 
Recall that under Qx, dBf = dBf + l3lxt=xdt. Now we can compute the expectation 
of Rq: 

Qx [Ro] = (3' / V(^s,Xs)(u,x^)du ds (25) 

JO JO 



To study this term, let us define the quantity 



'0 Jo ny/\ogn{\s -u\ + l) 
and the random variable (under P) Y : 



du ds 



(26) 



Y :-- 



'■{\X,-Xu\<C4y/\s-u\} 

n^/logn{\s -u\ + l] 



du ds 



(27) 
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Roughly speaking, Dn is what Y would be if C4 = 00. It is easy to check, with the 
central limit theorem, that we can take C4 large enough such that P[Y] > (1 — 6)Dn, 
uniformly in n. 

Then, as F < Dn a.s. it implies that 



P. 



Y < 



< 26 



(28) 



and a quick computation of the value of Dn gives us : 



Dn > v^logn. 



If we inject this result in the term in (|25|) , we have that with large P^j-probability, 
Now we can use the only free hand we have left : the value of n 



(29) 



let n = exp 



then Qx [Ro] > 



/a 



> 2 exp 



(30) 



if we take C5 large enough. 

Now we have to bound the variance of Rq under Qx- We can decompose Rq, 



„^j„Jo \Jo J Jo Jo 



{s,Xs){u,Xu) 



du ds 



x,yeJo 

Jo 



JO 



V^s,x){u,x^)du dB^^ + J2 



n / rs 



\J0 



{s,x){u,y)dBl ) dB'l 



Now, using that {x + y Y < 2a;^ + 2?/^, we can bound the variance by 

(l"n \ 2 pn PS 

s,x){u,Xu)du 1 + 2 2^ / / ,^)C/m rfs. 

-'0 / ,„^,.,,,^tJo Jo 



x£Jo ye Jo ' 



We already computed the right-hand term when we studied V^arg (_Ro)/ it is < 1. With 
similar techniques we find the other one : 



El Vis^X^x.^du < j^^^, 

sup / V(s,x){u,Xu)du < 

x&jn Jo 



therefore 2l3^n ^ ( j V(^s,x){u,x^)du J < 
IS Jo 



2/32 



100C| 



< 1 
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Now, using (|24l) , (|28|) and (|30|) , we can compute (|23|) 



exp(/;^(i?o)) < Po 



max \Xs\ > (C3 — l)\/n 

s&[0,n] 



P. 



(31) 



<25 + e-^' + 3e-2^' 

so that our result is proved provided that K and 5 have been chosen respectively large 
and small enough. 



In conclusion, we have proved (|23l) , which with (|TT|) , (|2T|) and (|22|) allow us to say 

and with ([l]) and the value of n from (|30)) : 

1 — m —2 /— c 

^g[log^T]<^ = - = exp (- 
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